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Abstract 

Large time behavior of solutions to abstract differential equations 
is studied. The results give sufficient condition for the solution to an 
abstract dynamical system (evolution problem) not to exibit chaotic 
behavior. The corresponding evolution problem is: 

u = A(t)u + F(t,u) + b(t), t > 0; u(0) = u . (*) 

Here u := u = u(t) £ H, H is a Hilbert space, i £ 1 + := [0, oo), 
A(t) is a linear dissipative operator: Ke(A(t)u, u) < —~f(t)(u,u), F(t,u) 
is a nonlinear operator, [|_F(f,ii)|| < co||it|| p , p > 1, co,p are positive 
constants, ||6(t)|| < (3(t), (5{t) > is a continuous function. 

Sufficient conditions are given for the solution u(t) to problem (*) 
to exist for all t > 0, to be bounded uniformly on M + , and a bound 
on is given. This bound implies the relation lim^oo = 

under suitable conditions on j(t) and f3(t). 

The basic technical tool in this work is the following nonlinear in- 
equality: 

g(t) < -l(t)g(t) + a(t, g{t)) + (3(t), t > 0; g(0) = g , 

which holds on any interval [0, T) on which g(t) > exists and has 
bounded derivative from the right, g{t) := lims^+o 9 ^ t+s } j — . It is 
assumed that 7(t), and /3(t) are real- valued, continuous functions of t, 
defined on R + := [0, oo), the function a(t,g) is defined for all t € M+, 
locally Lipschitz with respect to g uniformly with respect to t on any 
compact subsets [0, T], T < oo. 

If there exists a function fi(t) > 0, fi(t) G C 1 (R+), such that 

a { f -W)) +m -lk( l{t) -W))- V( - 0; " (0,9(0) ^ 

then g(t) exists on all of K+, that is T = oo, and the following estimate 
holds: 

< git) < Vt > 0. 

If /x(0)5(0) < 1, then < g{t) < Vt > 0. 
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1 Introduction 

A classical area of study is stability of solutions to evolution problems. We 
identify an evolution problem with an abstract dynamical system. An evolu- 
tion problem is described by an equation 

u(t) = F 1 (t,u), u(0)=u . (1) 

Here F\ : X — > X is a nonlinear operator in a Banach space X , u = u(t) = ^ . 
Quite often it is convenient to assume A to be a Hilbert space H, because 
the energy is often interpreted as a quantity (u, u) in a suitable Hilbert space. 
Suppose that F^t, 0) = and u = 0. Then u = is a solution to (JTJ). 
A. M. Lyapunov in 1892 published a classical work on stability of motion, 
where he studied equation (CQ) in the case X = MJ 1 and F\ analytic function 
of u. If Fx(t, 0) = 0, and F\ is twice Frechet differentiate, then one can 
write Fi(t,u) = A(t)u + F(t,u), where A(t) is a linear operator in X and 
\\F(t, u)\\ = 0(||w|| 2 ), ||w|| — >■ 0. This representation is a linearization of F 1 
around the point u = 0. Lyapunov defined the notion of stability (Lyapunov 
stability) of the equilibrium solution u = towards small perturbation of the 
data uq. He calls this solution stable (Lyapunov stable), if for any e > there 
is a 5 = 5(e) > such that if inequality ||m || < 5 holds then sup t>0 ||w(t)|| < e. 
Note that this definition implies the global existence of the solution to problem 
([T]) for all Uq in the ball ||m || < 5. 

The equilibrium solution u = is unstable if it is not Lyapunov stable. 
This means that there is an e > such that for any 5 > there is a uq, 
|| wo || < 5, and a t$ > such that ||ii(i,5)|| > e. 

One can give similar definitions for stability and instability of a solution to 
problem ([1]) with u ^ 0. In this case one calls the solution u — u(t; u ) stable 
if all the solutions u(t; w ) to problem ([1]), with w in place of u , exist for all 
t > and satisfy the inequality sup f>0 \\u(t;u ) — u(t;w )\\ < e provided that 
||wo — u>o|| < S. 

A solution u(t; uq) is called asymptotically stable if it is stable and there 
is a 5 > such that all the solutions u(t; wq) with \\u — Wo\\ < 8 satisfy the 
relation lim^oo \\u(t; uq) — u(t; Wo)\\ = 0. 

The equilibrium solution u = is asymptotically stable if it is stable and 
there is a 5 > such that all the solutions u(t; uq) with ||n || < 5 satisfy the 
relation lim^oo \\u(t; u )\\ = 0. 

Consider problem (JT]) with Fiit^u) + 4>(t,u) in place of Fi(t, u). The term 
(p(t, u) is called persistently acting perturbations. The equilibrium solution 
it = is called stable with respect to persistently acting perturbations if for 
any e > there exists a 5 = 5(e) > such that if \\(j)(t, u)\\ < 5 and ||wo|| < ^ ; 
then sup 4>0 \\u(t; u ) \\ < e. 

Stability of the solutions and their behavior as t — )■ oo are of interest in 
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a study of dynamical systems. For example, if the equilibrium solution is 
asymptotically stable, then it does not have chaotic behavior. 

If A(t) = A is independent of time and X = IR n , then Lyapunov obtained 
classical results on the stability of the equilibrium solution to problem ([1]). He 
assumed that F is analytic with respect tow £ M n , so that \F(t, u)\ < c\u\ 2 in a 
neighborhood of the origin, and c > is a constant. Lyapunov has proved that 
if the spectrum cr(A) of A lies in the half-plane Rez < 0, then the equilibrium 
solution u = is asymptotically stable, and if at least one eigenvalue of A lies 
in the half-plane Rez > 0, then the equilibrium solution is unstable. 

If some of the eigenvalues of A lie on the imaginary axis and F = 0, so that 
problem ([T]) is linear, and if the corresponding Jordan cells consist of just one 
element, then the equilibrium solution is stable. Otherwise it is unstable. 

Thus, a necessary and sufficient condition for Lyapunov stability of the 
equilibrium solution of the linear equation u = Au in M n is known: the spec- 
trum of A has to lie in the left complex half-plane: a C {Z : Rez < 0}, and 
the Jordan cells corresponding to purely imaginary eigenvalues of A have to 
consist of just one element. 

If F 0, then, in general, when the spectrum of A lies in the left half 
plane of the complex plane, and some eigenvalues of A lie on the imaginary 
axis, the stability cannot be decided by the linearized part A of F\ only. One 
can give examples of A such that the nonlinear part F can be chosen so that 
the equilibrium solution u = is stable, and F can also be chosen so that this 
solution is unstable. For instance, consider u = cu 3 , where c = const. This 
equation can be solved analytically by separation of variables. The result is 
u(t) = [u- 2 (0) - 2ct]-°- 5 . Therefore, if c < and |«(0)| < 5, 5 > 0, then the 
solution exists for all t > 0, and is asymptotically stable. But if c > 0, then the 
solution blows up at a finite time t&, the blow-up time, and tb = [2cit 2 (0)] _1 . 
In this case the zero solution is unstable. 

If A = A(t) the stability theory is more complicated. The case of periodic 
A(t) was studied much due to its importance in many applications (see |27J). 

The stability theory in infinite-dimensional spaces, for example, in Hilbert 
and Banach spaces, was developed in the second half of the 20-th century, see 
[Tj and references therein. Again, the location of the spectrum of A(t) plays 
an important role in this theory. 

The basic novel points of the theory presented below include sufficient con- 
ditions for the stability and asymptotic stability of the equilibrium solution 
to abstract evolution problem (JTJ) in a Hilbert space when o~(A(t)) may lie in 
the right half-plane for some or all t > 0, but sup o~(ReA(t)) — > as t — > oo. 
Therefore, our results are new even in the finite-dimensional spaces. 

The technical tool, on which our study is based, is a new nonlinear differ- 
ential inequality. The results are stated in several theorems and illustrated by 
several examples. These results are taken from the cited papers by the author 
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(see [14]- [22]), and, especially, from paper [20] ■ In the joint papers by the 
author's student N. S.Hoang and the author one can find various additional 
results on nonlinear inequalities (see [5]- [10]). Some versions of this inequality 
has been used in the monographs [12] and [13], where the Dynamical Systems 
Method (DSM) for solving operator equations was developed. 

The literature on stability of solutions to evolution problems and their 
behavior at large times is enormous, and we refer the reader only to the papers 
and books directly related to the novel points mentioned above. 

Consider an abstract nonlinear evolution problem 

it = A(t)u + F(t,u) + b(t), ^ := ^' ( 2 ) 

«(0) = u , (3) 

where u(t) is a function with values in a Hilbert space H, A(t) is a linear 
bounded dissipative operator in H, which satisfies inequality 

Re(A(t)u,u) < -7(t)||u|| 2 , t > 0; VueH, (4) 

where F(t,u) is a nonlinear map in H, 

\\F(t,u)\\<c \\u(t)r, p>l, (5) 

\\b(t)\\ <P(t), (6) 

j(t) > and j3(t) > are continuous function, defined on all of M + := [0, oo), 
c > and p > 1 are constants. 

Recall that a linear operator A in a Hilbert space is called dissipative if 
Ke(Au,u) < for all u e D(A), where D(A) is the domain of definition of A. 
Dissipative operators are important because they describe systems in which 
energy is dissipating, for example, due to friction or other physical reasons. 
Passive nonlinear networks can be described by equation (T5]) with a dissipative 
linear operator A(t), see [18], [19], Chapter 3, and [26] . 

Let a := a(A(t)) denote the spectrum of the linear operator A(t), IT := 
{z : Rez < 0}, £ := {z : Rez = 0}, and p(u,£) denote the distance between 
sets a and I. We assume that 

o C IT, (7) 

but we allow limt-Kx, p{a, £) = 0. This is the basic novel point in our theory. 
The usual assumption in stability theory (see, e.g., [1]) is sup zer7 Rez < —70, 
where 70 = const > 0. For example, if A(t) = A*(t), where A* is the adjoint 
operator, and if the spectrum of A(t) consists of eigenvalues \j(t), > \j(t) > 
\j + i(t), then, we allow lim^ooA^t) = 0. This is in contrast with the usual 
theory, where the assumption is Ai(t) < —70, 70 > is a constant, is used. 
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Moreover, our results cover the case, apparently not considered earlier in 
the literature, when Ke(A(t)u,u) < -y(t) with ^(t) > 0, lim t _ >00 7(t) = 0. This 
means that the spectrum of A(t) may be located in the half-plane Rez < j(t), 
where j(t) > 0, but lim t ^ 00 7(t) = 0. 

Our goal is to give sufficient conditions for the existence and uniqueness of 
the solution to problem ©-([3]) for all t > 0, that is, for global existence of u(t), 
for boundedness of sup i>0 < oo, or to the relation lim^oo ||w(t)|| = 0. 

If b{t) = in (|2D, then u(t) = solves equation Q2J and u(0) = 0. This 
equation is called zero solution to @ with b(t) = 0. 

Recall that the zero solution to equation (j2J) with b(t) = is called Lya- 
punov stable if for any e > 0, however small, one can find a 5 = 5(e) > 0, 
such that if ||w || < 5, then the solution to Cauchy problem (J2])-([2]) satisfies 
the estimate sup t>0 ||w(t)|| < e. If, in addition, lim^oo = 0, then the 

zero solution to equation ([7]) is called asymptotically stable. 

If b(t) ^ 0, then one says that ©-© is the problem with persistently 
acting perturbations. The zero solution is called Lyapunov stable for problem 
(EJ)- (EJ) with persistently acting perturbations if for any e > 0, however small, 
one can find a 5 = 5(e) > 0, such that if ||tt || < 5, and sup t>0 \\b(t)\\ < 5, then 
the solution to Cauchy problem (J2])-© satisfies the estimate sup t>0 ||w(£)|| < e. 

We do not discuss here the method of Lyapunov functions for a study of 
stability (see, for example, [11], [25]). 

The approach, developed in this work, consists of reducing the stability 
problems to some nonlinear differential inequality and estimating the solutions 
to this inequality. 

In Section 2 the formulation and a proof of two theorems, containing the 
result concerning this inequality and its discrete analog, are given. In Section 3 
some results concerning Lyapunov stability of zero solution to equation ([2]) are 
obtained. In Section 4 we derive stability results in the case when 7(i) > 0. 
This means that the linear operator A(t) in (j2J) may have spectrum in the 
half-plane Re^ > 0. 

Our results are closely related to the Dynamical Systems Method (DSM), 
see [2], [7], [8], [To]. Recently these results were applied to biological problems 
Q23J) and to evolution equations with delay (|22j). 

In the theory of chaos one of the reasons for the chaotic behavior of a 
solution to an evolution problem to appear is the lack of stability of solutions 
to this problem ([2], [3]). The results presented in Section 3 can be considered 
as sufficient conditions for chaotic behavior not to appear in the evolution 
system described by problem (J2])-(J3J). 
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2 A differential inequality 

In this Section an essentially self-contained proof is given of an estimate for 
non-negative solutions of a nonlinear inequality 

git) < -l(t)g(t) + a(t, git)) + 0{t), t > 0; g(0) = g ; 9'=^- (8) 

In Section 3 some of the many possible applications of this estimate (see esti- 
mate ffT2l) below) are demonstrated. 

It is not assumed a priori that solutions git) > to inequality (jSJ) are 
defined on all of M + , that is, that these solutions exist globally. In Theorem 1 
we give sufficient conditions for the global existence of g(t). Moreover, under 
these conditions a bound on g(t) is given, see estimate (I12p in Theorem 1. This 
bound yields the relation lim^oo git) = if lim^oo ji{t) = oo in (fl2l) . 

Let us formulate our assumptions. We assume that git) > 0. We do not 
assume that the functions 7, a and /3 are non-negative. However, in many 
applications the functions a and /3 are bounds on some norms, and then these 
functions are non- negative. The function j(t) is often (but not always) non- 
negative. For example, this happens if ^(t) comes from an estimate of the type 
(Au,u) > 7(u,u). If the functions a and f3 are bounds from above on some 
norms, then one may assume without loss of generality that these functions are 
smooth, because one can approximate a non-smooth function with an arbitrary 
accuracy by an infinitely smooth function, and choose this smooth function to 
be greater than the function it approximates. 

Assumption A x ). We assume that the function git) > is defined on some 
interval [0,T), has a bounded derivative g(t) := lim s ^+o 9 ^ +5 ^~ g ^ from the 
right at any point of this interval, and git) satisfies inequality (IE]) at all t at 
which git) is defined. The functions j(t), and (3(t), are real- valued, defined on 
all of M + and continuous there. The function a(t, g) is continuous on M + x M + 
and locally Lipschitz with respect to g. This means that 

\a{t,g)-a{t,h)\<L{T,M)\g-h\, (9) 

if t G [0, T], \g\ < M and \h\ < M. Here M = const > and L(T, M) > is a 
constant independent of g, h, and t. 

Assumption A2). There exists a C 1 (M + ) function fi(t) > 0, such that 

4^) +w -^( 7(t) -ID' Vi -°- (io) 

and 

M%(o) < 1. (11) 

One can replace the initial point t = by some point to G and assume 
that the interval of time is [to,^o + T), and that inequalities hold for t > to? 
rather than for t > 0. The proofs and the conclusions remain unchanged. 
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Theorem 1. If Assumptions A i) and A%) hold, then any solution g(t) > 
to inequality (jBJ) exists on all ofWL + , i.e., T = oo, and satisfies the following 
estimate: 

0<g(t)<-^ VteR+. (12) 

If 12(0)9(0) < 1, then < g(t) < ^ V£ G R + . 
Remark 1. 7/ lim^oo fi(t) = oo, then lim^oo g(t) = 0. 
Proof of Theorem [IJ Let us rewrite inequality for \i as follows: 

- i{t)»-\t) + a(t, fi-\t)) + 0(t) < (13) 
Let 0(t) solve the following Cauchy problem: 

0(t) = - 7 (t)0(t) + a(t,0(t)) + /3(t), t>0; 0(0) = O . (14) 

The assumption that a(t,g) is locally Lipschitz with respect to g guarantees 
local existence and uniqueness of the solution 0(t) to problem (fl4j) . From the 
comparison result (see A Comparison Lemma proved below) it follows that 

<P(t)<^\t) Vt>0, (15) 

provided that 0(0) < /i _1 (0), where <p(t) is the unique solution to problem 
(TT51) . Let us take 0(0) = g(0). Then 0(0) < /i _1 (0) by the assumption, and an 
inequality, similar to ffT5|) . implies that 



g(t)<<p(t) te[o,T). (16) 

Inequalities 0(0) < /i -1 (0), f|T5|) . and ffTBT) imply 

^)<0(t)<//- 1 (t), tG[0,T). (17) 

By the assumption, the function /i(t) is defined for all t > and is bounded 
on any compact subinterval of the set [0, oo). Consequently, the functions 0(t) 
and g(t) > are defined for all t > 0, and estimate ffl2|) is established. 

If g(0) < /i~ 1 (0), then one obtains by a similar argument the strict inequal- 
ity g{t) </i" 1 (t), t>0. 

Theorem [1] is proved. □ 
Let us now prove the comparison result that was used above, see, for ex- 
ample, jl], Theorem III. 4.1. 
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A Comparison Lemma. Let 

<j>(t) = f(t,<f>), 0(0) = 0o, (*) 

and 

ij)(t)=g(t,il)), ^(0) =Vo- (**) 

Assume ip > O , and 

g(t,x)> f(t,x) (***) 

/or any t and x /or which both f and g are defined. Assume that f and g 
are continuous functions in a set [0, s) x (a, 6), 0o € (a, 6), -0 zs t/ie maximal 
solution to (**) and is any solution to (*). Then (f>(t) < ip{t) on the maximal 
interval [0, T) of the existence of both and ip. 

Proof of the Comparison Lemma. First, let us assume for simplicity that 
problems (*) and (**) have a unique solution. Later we will discard this 
simplifying assumption. If / and g satisfy a local Lipschitz condition with 
respect to 0, respectively, ip, then our simplifying assumption holds. Assume 
secondly, also for simplicity, that g(t,x) > f(t,x). Under this simplifying 
assumption it is easy to prove the conclusion of the Lemma, because the graph 
of -0 must lie above the graph of -0 for t > 0. Indeed, in a small neighborhood 
[0, S), where 5 > is sufficiently small, the graph of ip lies above the graph of 
0. This is obviously true if O < -0o> because of the continuity of and -0. If 
0o = ipo, then the graph of -0 lies above the graph of because 0(0) < -0(0) 
due to the assumption /(O,0 O ) < g(O,0 o ) = g(O,-0o)- To check the last claim 
assume that there is a point t\ G [0, T) such that 0(ti) = "0(£i), and (pit) < ip(t) 
for t e (0,ti). Then 0(t) - 0(i x ) < ip(t) - 0(ti). Divide this inequality by 
t — ti < and get 

0(t)-0(tQ ^(t)-^( tl ) 

t-h t-t± 

Pass to the limit t — > t 1: t < t 1: in the above inequality, use the differential 
equations for and ip and the equality 0(ti) = and obtain the following 

relation: 

f(t U 0(t X )) = 0(ti) > #0 = V(ii)) = ff(tl, 0(tl))- 

This relation contradicts the assumption f(t,x) < g{t,x). The contradiction 
proves the conclusion of the Comparison Lemma under the additional assump- 
tion f(t,x) < g(t,x). 

To prove the Comparison Lemma under the original assumption f(t,x) < 
g(t, x), let us consider problem (*) with / replaced by /«:=/ — - < /• Let n 
solve problem (*) with / replaced by /„, and with the same initial condition 
as in (*). Since f n (t, x) < g(t, x), then, by what we have just proved, it follows 
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that 4> n (t) < if>(t) on the common interval [0, T n ) of the existence of <fr n and 
if). By the standard result about continuous dependence of the solution to (*) 
on a parameter, one concludes that lim^ooT^ = T and lim^oo 4> n (t) = <p(t) 
for any t G [0, T). Therefore, passing to the limit n — > oo in the inequality 
4>n{t) < ip{t) one gets the conclusion of the Comparison Lemma under the 
original assumption f(t,x) < g(t,x). 

If the simplifying assumption concerning uniqueness of the solutions to (*) 
and (**) is dropped, then (*) and (**) may have many solutions. The limit of 
the solution <p n is the minimal solution to (*). If one considers problem (**) 
with g replaced by g n := g + - > g, and denotes by ip n the corresponding 
solution, then the limit \\va. n ^ ) . oa 'ip n {t) = ip(t) is the maximal solution to (**). 
In this case the above argument yields the conclusion of the Lemma with ip(t) 
being the maximal solution to (**), and <fi(t) being any solution to (*). The 
Comparison Lemma is proved. □ 

Remark 2. // <p(t) is bounded from below for all t > 0, so that c < <p(t) 
for all t > 0, and if)(t) exists globally, that is, for all t > 0, then the inequality 
c < < ip{t) and the continuity of f(t,x) on the set [0, oo) x IR imply that 
any solution <fi to (*) exists globally. Indeed, if it would exist only on a finite 
interval [0, T) then it has to tend to infinity as t — > T, but this is impossible 
because the bound c < (f>(t) < ip{t) and the global existence and continuity of 'if) 
do not allow <p(t) to grow to infinity as t — » T. 

Let us formulate and prove a discrete version of Theorem [TJ 

Theorem 2. Assume that g n > 0, a(n,g n ) > 0, 

g n +i < (1 - h n j n )g n + h n a(n, g n ) + h n (3 n ; h n >0, < h n ^ n < 1, (18) 
and a(n,g n ) > a(n,p n ) if g n > p n . If there exists a sequence fi n > such that 

a(n, -)+Pn< -(Tn - ^l 1 "^ ), (19) 

f^n fJ'n >l"n 

and 

1 

go < —, 
Ho 

then 

o<g n <—, Vn > 0. 

/in 

Proof. For n = inequality (12"T|) holds because of (|2U|) 
for all n < m and let us check that then it holds for n = 
Theorem [2] is proved. 

Using the inductive assumption, one gets: 

g m +i < (1 - ^m7m)— + h m a(m, — ) + h m (3 m . 

/*m 



(20) 
(21) 

. Assume that it holds 
= m + 1. If this is done, 

□ 
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This and inequality ( 1191) imply 



< ('\ h \ -I- h ^ ( — ^' m + 1 ^ r 



i •"■rrt v i in t 



The last inequality is obvious since it can be written as 

-(/im - /Wl) 2 < 0. 

Theorem [2] is proved. □ 

Theorem [5] was formulated in [3] and proved in [0] . We included for com- 
pleteness a proof, which is shorter than the one in [BJ. 

Let us give a few simple examples of applications of Theorem 1. 
Example 1. Consider the inequality 

g(t)<tg~(t + l) 2 g 2 -2(t + l)- 2 . (22) 

Assume g > 0. Choose fi(t) = t + 1. Then inequality fllUp holds if 

(t + l)[-(t + l) 2 (i + I)" 2 - 2(t + I)" 2 ] <-t-(t+ l)" 1 , 

and (?(0) < 1. Thus, inequality ( fTUj) holds if 

-t - 1 - 2(t + l)" 1 < -t - (t + 

This inequality holds obviously. Therefore, any g > 0, that satifies inequalities 
( 122|) and (7(0) < 1, exists for alH > and satisfies the estimate 

< g{t) < ' 



t+l 

In this example the linearized problem 

g(t) =tg-2(t + l)- 2 , g{0)=g , 
has a unique solution 



g{t) = e*' 2 \g{0)-2 f c "4 ( s + 1)-^]. 



This solution tends to infinity as t — > 00. 
Example 2. Consider a classical problem 



= A(t)u + F(t,u), u(0) = u , (23) 



12 



where A(t) is a linear operator in R n and F is a nonlinear operator. Assume 
that Ke(A(t)u, u) < —j(u,u), where 7 = const > 0, and \\F(t, u)\\ < c\\u\\ p , 
p = const > 1, c = const > 0, and 1 1 ■ 1 1 is the norm of a vector in M n . We also 
assume that equation ([23]) has the following property: 

Property P: If a solution to ( 123]) is defined on the maximal interval of its 
existence [0,T) and T < 00, then lim^j^o = 00. 

It is known (see, for example, that Property P holds if F(t,u) is a 
continuous function on [0, T] x M n . 

By Peano's theorem the Cauchy problem 

ii(t) = f(t,u), u(0)=u o , (24) 

where u G M n , has a local solution on an interval [0,a), provided that / is 
a continuous function on [0,T] x D(uo), where a G (0,T) and -D(«o) is a 
neighborhood of uq. This solution is non-unique, in general. One can give 
an explicit estimate of the length a of the interval on which the solution does 
exist. Namely, a = min(T, -p), where M := max\ u - Uo \<b t te[o,T]\f{t, u)\, and the 
neighborhood D{u$) is taken to be the set {u : \u — u \ < b}. 

It is known that in every infinite- dimensional Banach space the Peano 
theorem fails. Therefore, in an infinite-dimensional Banach space we assume 
that problems f l2"3"]) and f ]2~4"]) have a solution, and if [0,T) is the maximal 
interval of the existence of the solution, then Property P holds. This happens, 
for example, if f(t, u) satisfies a local Lipschitz condition with respect to u and 
is continuous with respect to t G [0,T]. Indeed, if a local Lipschitz condition 
holds, then the local interval of the existence of the solution to the Cauchy 
problem ( )24j) is of the length b = min(RM~ l , L), provided that / is continuous 
with respect to t and satisfies the estimates | |/(t, u)\ \ < M, \\f(t, u)— f(t, v)\\ < 
L\ \u — v\ |, in the region [0, T] x B(uo, R), B(uo, R) := {u : \\u — Uo\\ < R}- 
Under these assumptions the solution to problem (124"]) is unique and stays in 
the ball B(u Q ,R) for t G [0,6]. 

To see that Property P holds for problem (12"41) if / satisfies a local Lipschitz 
condition with respect to u, assume that the solution to (T241) does not exist 
for t > T. Under our assumptions, if the solution u of problem (]24]) satisfies 
the inequality sup 0<i<T ||w(t)|| < 00, then the constants M, L and R are finite. 
Therefore b > 0. Take the initial point to = T — 0.5b. By the local existence 
theorem the solution u(t) exists on the interval [T — 0.5b, T + 0.56]. This is 
a contradiction, since we have assumed that this solution does not exist for 
t > T. This contradiction proves that Property P holds for problem (1241) if / 
satisfies a local Lipschitz condition. 

Let us use Theorem 1 to prove asymptotic stability of the zero solution 
to (125]) and to illustrate the application of our general method for a study 
of stability of solutions to abstract evolution problems, the method that we 
develop below. 
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Let g(t) := ||tt(t)||, where the norm is taken in R n . Take a dot product 
of equation (1231) with u, then take the real part of both sides of the resulting 
equation and get 

Re(u, u)= gg = Re(Au, u) + Re(F(t, u),u) < -jg 2 + cg p+1 . 

Since g > 0, one obtains from the above inequality an inequality of the type 
(jHJ), namely, 



where A = const > 0, a = const G (0, 7). Note that a can be chosen arbitrarily 
close to 7. We choose A later. Denote b := 7 — a > 0. Then inequality f fTTj) 
holds for any g(0) if c > is sufficiently small. Inequality f lTUj) holds if 



Since p > 1 this inequality holds if cA _( - p_1 ^ < 6. In turn, the last inequality 
holds for an arbitrary fixed c > and an arbitrary small fixed b > provided 
that A > is sufficiently large. 

One concludes that for any initial data Uq the solution to f[2"3"j) exists globally 
and admits an estimate \ \u(t)\\ < A _1 e~ a ' ; where the positive constant a < 7 
can be chosen arbitrarily close to 7 if the positive constant c is sufficiently 
small. 

The above argument remains valid also for unbounded, closed, densely de- 
fined linear operators A(t), provided that Property P holds. 

If A(t) is a generator of a Co semigroup T(t), and F satisfies a local 
Lipschitz condition, then problem ff23l is equivalent to the equation u = 
T(t)F(t,u), and this equation may be useful for a study of the global exis- 
tence of the solution to problem (I23"|) (see [24]). 

Example 3. Consider an example in which the solution blows up in a 
finite time, so it does not exist globally. Consider the problem 

ii-Au = u 2 in [O.oojxDci™; u N = 0; u(0,x) = u {x). (25) 

Here D is a bounded domain with a smooth boundary S, N is an outer unit 
normal to S, uq > is a smooth function. Let 



Integrate equation (|2"5"|) over D and get g(t) = f D u 2 dx. Use the inequality 




p = const > 1, 



< 7 — a = b. 
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where c = c(D) = const > 0, and get g > g 2 /c. Integrating this inequality, 
one obtains g(t) > — ct] _1 . Since c > and g > it follows that 

lim <7(t) = oo, 

where £& := — is the blow-up time, and t < t b . Consequently, for any initial 
data with g > the solution to (j25p does not exist globally. 
Example 4. Consider the following equation 

u + A[t)u + c/)(u) - = fit, u), u(0, x) = u (x), (26) 

where u = u(t,x), <fi and ip(t,u) are smooth functions growing to infinity as 
\u\ — > oo. Let us assume that 

u<p{u)>0, mp{t,u)>0 Vt>0, 

and 

uip(t,u) < a(t)\u\ 3 , \uf(t,u)\ < P(t)\u\, 

where a(t) > and (3{t) > are continuous functions, x G D C M n , D is a 
bounded domain, 

Re(Au,w) > 7(m,m) Vm G -D(^4), 7 = const > 0, 

A is an operator in a Hilbert space = L 2 (D), the domain of definition of A, 
D(A), is a dense in if linear set, (u, v) is an inner product in H, | \u\ | 2 = (w, w). 
An example of A is A = —A, the Laplacean with the Dirichlet boundary 
condition on S, the boundary of D. Denote g(t) := We want to 

estimate the large time behavior of the solution u to fl26|) . 

Take the inner product in H of fl2"61) and u, then take real part of both 
sides of the resulting equation and get 

99 < -19 2 + ag 3 + fig. 

Since g > one obtains an inequality of the type ([H]), namely 

9 < -19 + a(t)g 2 + 0. 

Now it is possible to use Theorem 1. 

Choose a{t) = Xe kt , where A and k are positive constants, k < 7. Assume 
that Ago < 1) where g := ||m (x)||. Then inequality ffTTl) holds for any initial 
data Mo, that is, for any g , if A is sufficiently small. Inequality (jTOl) holds if 

a{t) f kt + Xe kt (3(t) < 7 - k. 
A 
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One can easily impose various conditions on a and (3 so that the above in- 
equality hold. For example, assume that a decays monotonically as t grows, 
< (7 — fc)/2, and (3{t) < z/e~ fc '*, where k' > k, k' — const, v > is a 
constant, \u < (7 — k)/2. Then inequality ffTUj) holds, and it implies that 



so that the exponential decay of as t — > 00 is established. 

In Sections 3 and 4 some stability results for abstract evolution problems 
are presented in detail. These results are formulated in four theorems. The 
basic ideas are similar to the ones discussed in examples in this Section, but 
new assumptions and new technical tools are used. 

3 Stability results 

In this Section we develop a method for a study of stability of solutions to 
the evolution problems described by the Cauchy problem (J2])- (J3J) for abstract 
differential equations with a dissipative bounded linear operator A(t) and a 
nonlinearity F(t,u) satisfying inequality ([5]). Condition ([5]) means that for 
sufficiently small the nonlinearity is of the higher order of smallness 

than ||-u(t)||. We also study the large time behavior of the solution to problem 
(EJ)-([3|) with persistently acting perturbations b{t). 

In this paper we assume that A(t) is a bounded linear dissipative operator, 
but our methods are valid also for unbounded linear dissipative operators A(t), 
for which one can prove global existence of the solution to problem (T5])-([3]). We 
do not go into further detail in this paper. 

Let us formulate the first stability result. 

Theorem 3. Assume that Re(Au,u) < — A;||-u|| 2 Vw G H , k = const > 0, 
and inequality (jlj) holds with ^(t) = k. Then the solution to problem (J2])- 
([3]) with b(t) = satisfies an esimate \\u(t)\\ = 0(e _ ^~^') as t — >■ 00. Here 
< e < k can be chosen arbitrarily small if ||wo|| is sufficiently small. 

This theorem implies asymptotic stability in the sense of Lyapunov of the 
zero solution to equation (T5]) with b(t) = 0. Our proof of Theorem 3 is new 
and very short. 

Proof of Theorem 3. Multiply equation ([2]) (in which b(t) = is assumed) 
by u, denote g = g(t) := take the real part, and use assumption (j4]) 

with j(t) = k > 0, to get 



—kt 



U(t)\\< 



A 



gg < -kg 2 + c g p+1 , 
If g(t) > then the derivative g does exist, and 



p > 1. 



(27) 
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as one can check. If g(t) — on an open subset of K+, then the derivative g 
does exist on this subset and g{t) = on this subset. If g(t) — but in in any 
neighborhood (t — 5, t + 5) there are points at which g does not vanish, then 
by g we understand the derivative from the right, that is, 

m := lim + ■)-»(') = lim ?*±i). 

This limit does exist and is equal to ||it(t)||. Indeed, the function u(t) is 
continuously differentiable, so 

lim 11"" + ^ , im MW±gWjl = | W) ||, 

s s 

The assumption about the existence of the bounded derivative g(t) from the 
right in Theorem 3 was made because the function does not have, in 

general, the derivative in the usual sense at the points t at which = 0, 

no matter how smooth the function u(t) is at the point r. Indeed, 

lim Kt+i)j| = Um MW+£Wj! = _ h 

s->-0 S s->— S 

because lim s ^_ ^ = — 1. Consequently, the right and left derivatives of ||n( 
at the point t at which ||w(i)|| = do exist, but are different. Therefore, the 
derivative of at the point t at which = does not exist in the 

usual sense. 

However, as we have proved above, the derivative g(t) from the right does 
exist always, provided that u(t) is continuously differentiable at the point t. 

Since g > 0, inequality ( 12T|) yields inequality ([8]) with j(t) = k > 0, 
f3(t) = 0, and a{t,g) = c g p , p > 1. Inequality (jTUI) takes the form 

-4 T <-Lffc-^y Vt>0. (28) 

Let 

H(t) = Xe bt , X,b = const > 0. (29) 

We choose the constants A and b later. Inequality ( |T0|) . with /x defined in ( 129~1) . 
takes the form 

A ,-i e C ( P -i) fcf + & < fc > V ^>°- (3°) 
This inequality holds if it holds at t = 0, that is, if 



(31) 
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Let e > be arbitrary small number. Choose b = k — e > 0. Then ( 13TT) holds 
if 

A>(-)^. (32) 
e 



Condition flU)) holds if 

IKII =0(0) < ~. (33) 

We choose A and b so that inequalities (l3~2i) and (!33l hold. This is always 
possible if b < k and ||«o|| is sufficiently small. 

By Theorem 1, if inequalities (l3~TT)-(l33l) hold, then one gets estimate ( TT2|) : 



0<g(t) = \\u(t)\\<^j^, Vt>0. (34) 

Theorem 3 is proved. □ 
Remark 3. One can formulate the result differently. Namely, choose 
A = ||«o|| _1 . Then inequality (J33l) holds, and becomes an equality. Substitute 
this A into ( 1311 and get 

colMI^ 1 + b<k. 

Since the choice of the constant b > is at our disposal, this inequality can 
always be satisfied if Cq\\uq\\ p ^ 1 < k. Therefore, condition 

c o\\ u o\\ p l < k 
is a sufficient condition for the estimate 

IK*) II < NI|e~ (fe ~ co||uo|rl)t , 



to hold (assuming that co||mo|| p_1 < k). 

Let us formulate the second stability result. 
Theorem 4. Assume that inequalities Qj-Qj hold and 

l(t) = ^ , qi < 1; ci, q x = const > 0. (35) 
Suppose that e G (0, Ci) is an arbitrary small fixed number, 
A>(^) 1/(P_1) and K0)||<i. 



TTien i/ie unique solution to ©-© ™^ K*) = exzsfo on a// o/R + and 

0<l|w(*)ll<w ^ , Vt>0. (36) 
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Theorem 4 gives the size of the initial data, namely, ||u(0)|| < v, for which 
estimate ( 1361) holds. For a fixed nonlinearity F(t,u), that is, for a fixed con- 
stant Co from assumption (jSJ), the maximal size of ||u(0)|| is determined by the 
minimal size of A. 

The minimal size of A is determined by the inequality A > (sa) 1 ^ ^ that 

is, by the maximal size of e G (0, C\). If e < c\ and ci — e is very small, then 
/ \i/(p-i) 

A > Xmin '■= ( ^ ) and A can be chosen very close to A min . 

Proof of Theorem 4- Let 

//(*) = A(l + *)", A, i/ = const > 0. (37) 

We will choose the constants A and v later. Inequality f JTOj) (with (3(t) = 0) 
holds if 

+ + TTt - (irV' W -°' (38) 

if 

gi < 1 and (p — l)v > gi, (39) 
then inequality ( l38l) holds if 

C ^T + " < ( 4 °) 



Ap- 

Let e > be an arbitrary small number. Choose 

u = ci - e. (41) 



Then inequality ([40]) holds if inequality (j32j) holds. Inequality fTTTj) holds be- 
cause we have assumed in Theorem 4 that ||w(0)|| < Combining inequalities 
( )32|) . fl33|) and (TT2]) . one obtains the desired estimate: 

0<\\u(t)\\=g(t)< ^ , Vt>0. (42) 

Condition ( 132|) holds for any fixed small e > if A is sufficiently large. Condi- 
tion ([33]) holds for any fixed large A if ||mo|| is sufficiently small. 

Theorem 4 is proved. □ 
Let us formulate a stability result in which we assume that b(t) ^ 0. The 
function b(t) has physical meaning of persistently acting perturbations. 
Theorem 5. Let b(t) ^ 0, conditions (j3])- (jSJ) and f[35l) ZioZd, and 

< - ° 2 , , (43) 
v y ~ (1 + t) 92 

where ci > and g 2 > are constants. Assume that 

q x < min{l, g 2 - z/, u(p - 1)}, ||u(0) || < A„ \ (44) 
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where Ao > is a constant defined in (ToTTl . and 

cl~ kp cl(p-l)r^— + v< Cl . (45) 
p — 1 

Then problem ©-(ED has a unique global solution u(t), and the following esti- 
mate holds: 

Proof of Theorem 5. Let g(t) := \\u(t)\\. As in the proof of Theorem 4, 
multiply fl2J) by w, take the real part, use the assumptions of Theorem 5, and 
get the inequality: 

g<- „ °\ g + c g p + - ° 2 . . (47) 

y ~ (l+t)9i y Uy (1 + t)«a V 1 

Choose fi(t) by formula (137)1 . Apply Theorem 1 to inequality ( 1471) . Condition 
(ITU)) takes now the form 

Co Aco ^ Ci , , . , , 

i — Ti — w T\ + t r + < t Vt > 0. 48 

If assumption (14*4)1 holds, then inequality ([4*8)1 holds provided that it holds for 
t = 0, that is, provided that 

C ° T + Ac 2 + z/< ci. (49) 



A/ J 

Condition O holds if 



0(0) < (50) 

The function h(X) := -jprr + Ac 2 attains its global minimum in the interval 
[0, oo) at the value 

A ^f(-iM 1/P , (B1) 



c 2 



and this minimum is equal to 



1 / \- P 

hmin = C%C 2 P (p-l)P 



p — 1 



Thus, substituting A = A in formula (|49)1 . one concludes that inequality (149)1 
holds if the following inequality holds: 

4^(p-l)?^- + i/< Cl , (52) 
p — 1 
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while inequality (TSUI) holds if 



||«(0)||<i (53) 
Therefore, by Theorem 1, if conditions ( 152|) - (|53|) hold, then estimate ( !T2|) yields 

IW()ll£ VTTF' V '-°' (54) 

where A is defined in (I5T1) . 

Theorem 5 is proved. □ 



4 Stability results under non-classical assump- 
tions 

Let us assume that Re(A(t)u,u) < 7(t)||u|| 2 , where j(t) > 0. This corresponds 
to the case when the linear operator A(t) may have spectrum in the right half- 
plane Re^ > 0. Our goal is to derive under this assumption sufficient conditions 
on j(t), a(t,g), and 0(t), under which the solution to problem (T5]) is bounded 
as t — > oo, and stable. We want to demonstrate new methodology, based on 
Theorem 1. By this reason we restrict ourselves to a derivation of the simplest 
results under simplifying assumptions. However, our derivation illustrates the 
method applicable in many other problems. 
Our assumptions in this Section are: 

/5(t) = 0, 7(t) = Cl (i + t)- m \ a(t,g) = c 2 (l + t)- m2 gP, p > 1. 

Let us choose 

^(t) = d + \(l + t)- n . 

The constants Cj,nij, X,d,n, are assumed positive. 

We want to show that a suitable choice of these parameters allows one to 
check that basic inequality f lTOj) for fi is satisfied, and, therefore, to obtain in- 
equality ( IT2l) for g(t). This inequality allows one to derive global boundedness 
of the solution to (T5]), and the Lyapunov stability of the zero solution to (T5]) 
(with uq = 0). Note that under our assumptions fi < 0, lim^oo /i(t) = d. We 
choose A = d. Then (2d) _1 < /i _1 (t) < dr 1 for all t > 0. The basic inequality 
(ITOl) takes the form 

Cl (l + t)- m ^ +C2 (l+ty m2 [d + \(l + t)- n ]- p+1 <nA(l + t)-"- 1 [rf + A(l+t)-"]- 1 , 

(55) 

and 

g {d + \)<l. (56) 
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Since we have chosen A = d, condition ( 1561) is satisfied if 

d=(2g y 1 . (57) 

Choose n so that 

n + 1 < minjm!, m 2 }. (58) 

Then (1551) holds if 

ci + c 2 cT p+1 < nAcT 1 . (59) 

Inequality (1591) is satisfied if ci and c 2 are sufficiently small. Let us formulate 
our result, which folows from Theorem 1. 

Theorem 6. If inequalities (1591) and (1581) hold, then 

< git) < [d + \{1 + ty 11 ]- 1 < d-\ Vt>0. (60) 

Estimate (160]) proves global boundedness of the solution u(t), and implies 
Lyapunov stability of the zero solution to problem (T5]) with b(t) = and 
m = 0. 

Indeed, by the definition of Lyapunov stability of the zero solution, one 
should check that for an arbitrary small fixed e > estimate sup <>0 ||w(£)|| < e 
holds provided that ||w(0)|| is sufficiently small. Let ||w(0)|| = g = 5. Then 
estimate fl60l) yields sup t>0 ||tt(t)|| < d~ x , and floTl) implies sup t>0 ||M(t)|| < 25. 
So, e = 25, and the Lyapunov stability is proved. □ 
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